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Abstract 

This  paper  deals  with  network  flow  analysis  and  its  application  to  the  calculation  of  flow  distribution  within  fuel  cells.  The  generated  node 
and  loop  equations  for  the  network  analysis  are  presented.  In  order  to  accommodate  changes  in  cell  design,  i.e.  changes  in  flow  network 
topology,  the  analysis  includes  automatic  loop  equation  generation  based  on  topology  analysis  and  network  search  algorithm.  The  calculated 
flow  distributions  were  in  excellent  agreement  with  experimental  data.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Flow  networks  exist  in  a  wide  range  of  practical  systems, 
including  transportation  systems,  water  distribution  sys¬ 
tems,  thermal  systems,  as  well  as  information  flow  through 
the  Internet.  In  a  thermal  system,  flow  network  analysis 
serves  a  critical  design  purpose  involving  fluid  flow,  heat  and 
mass  transfer  for  equipments  such  as  heat  exchanger,  boiler 
and  most  recently,  fuel  cell  stacks.  Analysis  of  the  flow 
distribution  provides  a  better  understanding  of  the  perfor¬ 
mance  of  such  network  flow  systems. 

In  a  fuel  cell  generating  system,  fuel  cells  are  connected 
electrically  in  series  to  create  a  fuel  cell  stack.  Usually,  fuel 
or  oxidant  flows  through  the  fuel  cell  in  multiple  channels 
arranged  in  a  complex  flow  network.  The  flow  distribution 
within  a  fuel  cell  stack  is  a  major  concern  facing  fuel  cell 
designers  since  it  has  a  significant  impact  on  fuel  cell 
performance  and  efficiency  [1],  Numerous  studies  dealing 
with  the  electrochemical  and  material  aspects  of  fuel  cells 
have  been  performed.  However,  to  reduce  manufacturing 
costs  and  improve  performance,  better  understanding  of  the 
fuel  and  oxidant  species  transport  processes  within  the  fuel 
cell  stack  is  critical  and  is  still  an  open  question.  Hirata  and 
Hori  [2]  have  studied  the  effects  of  the  gas  flow  uniformity  on 
cell  performance.  They  found  that  a  uniform  inlet  flow 
distribution  results  in  improved  performance  for  the  co-flow 
fuel  cell  configuration.  To  achieve  uniform  flow,  many 
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designs  often  use  a  complex  flow  network,  and  inevitably, 
increase  the  stack  cost,  (Sadowski  et  al.  [3]). 

Analysis  of  the  fuel  cell  stack  performance  requires 
knowledge  of  the  flow  distribution  within  the  stack.  For 
this  purpose,  computational  fluid  mechanics  (CFD)  presents 
a  helpful  tool.  This  method  is  limited  to  simple  geometric 
configurations  with  few  flow  channels.  However,  flow 
through  the  distribution  manifold  of  a  fuel  cell  generally 
forms  a  complex  flow  network  (Fig.  1),  which  makes  CFD 
computations  exceedingly  costly  and  time  consuming  due  to 
the  difficulty  of  generating  complex  computational  grid. 
Hence,  flow  network  analysis  offers  a  useful  tool  for  such 
applications.  Various  mathematical  techniques  and  compu¬ 
ter  software  have  been  developed  for  solving  the  problem  of 
flow  and  pressure  distributions  in  network  systems  [5], 
However,  these  methods  have  restrictions  and  cannot  cur¬ 
rently  be  used  to  estimate  the  fuel  cell  stack  performance. 
With  its  long  tested  history  in  solving  network  flow  analysis 
problems  for  thermal  systems,  SINDA/FLUINT  [6]  has  a 
great  potential  in  dealing  with  the  flow  distribution  issue 
within  fuel  cell  stack  designs.  SINDA/FLUINT  may  also 
be  extended  to  calculate  the  parameter  profiles  within  a  fuel 
cell  stack,  as  long  as  an  adequate  fuel  cell  model  can  be 
implemented  into  its  network  analysis  procedure. 

In  this  investigation,  methods  for  complex  network  flow 
analysis  have  been  developed  and  compared  with  experi¬ 
ments.  The  presentation  will  begin  with  simple  distribution 
systems  with  manually  generated  node  and  loop  equations. 
In  order  to  accommodate  design  changes,  i.e.  changes  in  the 
flow  network  topology,  an  automatic  scheme  for  loop  equa¬ 
tion  generation  based  on  topology  analysis  and  network 
search  algorithm  will  then  be  presented. 
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Nomenclature 

As  cross-section  area  (m2) 

Aw  wall  area  of  a  channel  (m2) 

b  the  number  of  branches 

Cf  frictional  coefficient 

d  the  number  of  nodes 

Du  hydraulic  diameter  (m) 

E  energy  term 

G  mass  flux  (kg/m2s) 

7/ser  the  energy  losses  through  the  loop 

//L  the  head  loss  in  the  branch 

K  the  loss  coefficient 

/  the  number  of  loops 

m  mass  flow  rate  (kg/s) 

For  p  pressure  (Pa) 

Re  Reynolds  number 

5  sign  function  of  flow  direction  with  respect 
to  node  (+ 1 ,  —  1  or  0) 

u  velocity  in  x  direction  (m/s) 

u  average  flow  velocity  in  a  channel  (m/s) 

x  x  coordinate  direction 

Subscripts 

i  branch  number 

j  node  or  loop  number 

Greek  symbols 

6  gas  channel  height  (m) 

p  dynamic  viscosity  of  solution  (N  s/m2) 

6  sign  function  for  branch  flow  direction 

(+1, -lorO) 

p  mass  density  (kg/m3) 

r  shear  stress  (N/m2) 

£,  pressure  loss  coefficient 


2.  Network  analysis 

A  complex  flow  system  can  be  reduced  into  a  graphical 
and  numerical  construction  as  illustrated  in  Fig.  2.  In  this 
graph,  nodes  and  branches  are  defined  and  numbered.  From 
the  graph  structure,  a  simple  loop  starting  from  one  node, 
passing  through  some  closed  branches,  and  returning  to  the 
original  node,  is  identified.  The  network  flow  distribution 
into  the  various  branches  can  be  computed  based  on  appli¬ 
cation  of  the  conservation  principles  at  the  nodes  and 
through  the  flow  loops. 

Consider  a  flow  network  as  shown  in  Fig.  1 .  The  funda¬ 
mental  equation  to  be  satisfied  is  mass  conservation  for  a 
node,  which  states  that  the  flow  into  and  out  of  each  node 
should  be  equal, 

y>,-  =  y (i) 

inlet  outlet 


or 


d+i 

ffSjiihi  =  0  (/=1,2  ,•••,</)  (2) 

i=  1 

where  ,S}(-  represents  the  sign  of  flow  in  branch  i  relative  to 
node  j,  and  the  subscript  i  represents  the  branch  number.  Sp  is 
equal  to  + 1  when  fluid  flows  into  the  node,  and  is  equal  to 
—  1  when  fluid  flows  out  from  the  node.  Sp  is  0  if  branch  i  has 
no  connection  with  node  j.  Here  d  is  the  total  number  of 
nodes,  while  /  is  the  total  number  of  loops. 

For  each  of  the  loops  in  the  network  (Fig.  1),  i.e.j  =  d  +  1, 
d  +  2,  ...,</  +  /,  the  energy  conservation  equation  can  be 
written  as 


HseI~J2E v  =  AE  O) 

where  HseT  is  the  energy  losses  through  the  loop,  Ep  is  the 
energy  input  to  the  fluids  by  pump  or  compressor,  and  A E  is 
the  difference  in  pressure  head  in  source  nodes. 

The  head  loss  in  each  series  channel  add  to  the  energy 
loss,  HseT,  that  is, 

Hser=  (4) 

where  //j  represents  the  head  loss  of  loop  j  resulting  from 
friction  associated  with  flow  through  branch  i.  It  can  be 
expressed  in  a  general  form  as 

Hr,jr  =  Y.  f  (5) 

where  Kt  is  the  loss  coefficient  for  different  head  losses  and 
the  exponent  n,  can  vary  among  the  components. 

For  the  interior  flow  distribution  within  a  fuel  cell  stack 
(see  for  example  Fig.  1),  no  energy  input  and  head  change 
through  the  node  is  present.  Hence,  the  energy  input  term, 
Ep,  and  the  pressure  difference  in  source  nodes,  A E,  can  be 
neglected.  Therefore,  from  Eq.  (3),  the  loop  or  energy 
equations  state  that  the  sum  of  pressure  changes  around 
any  loop  j  should  be  zero,  i.e. 

d+l 

ffdji  AP,  =  0  (j  =  d  +  1,  d  +  2, . . . ,  d  +  /)  (6) 

i=  1 

where  dp  is  a  sign  convention  representing  the  direction  of 
flow  in  branch  i  relative  to  loop  j.  dp  is  +1  when  fluid  flows 
in  the  same  direction  as  marked  in  the  figure  and  —  1  if  the 
fluid  flow  direction  is  opposite  to  the  marked  direction.  Qp  is 
0  if  branch  i  does  not  belong  to  loop  j. 

The  pressure  drop  of  branch  i  in  flow  loop  /  is  composed  of 
friction  loss  and  local  loss,  which  can  be  written  in  the  form 
of 

^4 

A  P,  =  AFiniet  +  A  P  outlet  +  APbend  +  t  —  (7) 

where  t  is  the  shear  stress,  which  can  be  found  from 


iij  A  v, 


(8) 


INLET  MANIFOLD 


SCHEMATIC  OF  CATHODE  USED  IN  GEORGIA  TECH  CODE 

Fig.  1.  Sample  diagram  for  a  complex  flow  network. 


Fig.  2.  Schematic  drawing  of  pressure  drop  experiment  set-up. 
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Cf  is  the  friction  coefficient,  which  is  dependent  on  the 
Reynolds  number  and  the  path  cross-sectional  geometry 
shape.  Here,  we  use  the  real  cross-sectional  area  As  and 
wall  area  Aw  to  account  for  the  flow  velocity  and  shear  force. 
Using  the  appropriate  friction  coefficient,  the  model  can  be 
tailored  to  any  flow  path  shape  and  thus  overcomes  the 
shortcomings  of  commercial  network  analysis  codes. 

The  transitional  loss  can  be  estimated  in  terms  of  the  loss 
coefficients 


P  —2 

A.Pinlet  =  inlet  ^  : 

P  —2 

A/V-nd  =  Cbcnd  ^ 


P  — 1 

APoutlet  =  t  outlet 


g  Ui  . 

'  s  inlet  2^  ^ 

(9) 

g  Ui  . 

-  Cbend 

GO) 

_  g  Ui  . 

Coutlet  2^  ^ 

(ID 

where  the  subscripts  inlet,  outlet  and  bend  represent  inlet, 
outlet  and  bend,  respectively.  The  loss  coefficients 
£iniet>  ^outlet!  and  £bend  for  various  types  of  junctions  are  listed 
in  Idelchik’s  handbook  [4].  The  transitional  loss  coefficient 
in  reference  [4]  is  expressed  as  a  function  of  the  connection 
geometry,  along  with  the  ratio  of  the  flow  rates  among 
branches.  It,  therefore,  accurately  includes  the  factors  affect¬ 
ing  the  local  pressure  loss,  which  is  frequently  assumed  to  be 
constant.  Substituting  Eqs.  (9)-(ll),  we  can  express  the 
pressure  drop  across  branch  i  in  the  following  form 


A  Pi  =  AP, 


inlet 


AP, 


outlet 


APbend  +  t  "T- 


(  A i n I c L  2^  +  ^outlet  2^  1  Cbend  ^ 2  ) 


(12) 


performed  for  generating  a  new  solution  until  the  difference 
between  the  new  solution  and  the  previous  solution  converge 
within  a  specified  tolerance.  The  global  linearization  method 
relaxes  the  solution  process  compared  to  a  solver  for  non¬ 
linear  system  of  equations.  The  network  analysis  program 
based  on  the  above  numerical  scheme  has  been  assigned  the 
name  of  “GT  code.” 

Following  the  procedure  described  above,  the  network 
equations  could  be  generated  manually.  However,  this 
would  limit  the  model’s  flexibility,  that  is,  if  the  flow 
network  were  to  be  modified  during  the  system  design 
process,  one  needs  to  redo  the  entire  process  manually. 
To  overcome  this  drawback,  the  evolution  program  method 
developed  by  Savic  and  Walters  [7],  was  adopted.  The 
hydraulic  solver  uses  a  method  based  on  loop  equations, 
namely,  the  linear  theory  method.  Network  topology  ana¬ 
lysis  is  performed  using  depth-first-search  and  breadth- 
first-search  algorithms  that  are  devised  to  ensure  the  con¬ 
nectedness  of  a  network  and  to  enable  easy  identification  of 
loops.  The  model  to  be  developed  will  enable  rapid  evalua¬ 
tion  of  a  network  for  different  topological  conditions  (i.e. 
the  connectivity  of  the  nodes  caused  by  different  piping 
configurations). 

The  automatic  network  analysis  code  written  in  Fortran90 
using  the  search  algorithm  [7]  has  been  tested.  The  auto¬ 
matically  generated  loop  equations  were  successfully  solved 
by  using  the  global  linearization  method.  The  results  from 
both  the  manually  generated  network  equation  and  the 
automatic  network  analysis  are  compared  with  experimental 
data,  as  well  as  the  result  from  SINDA/FLUINT. 


3.  Experimental  verification  of  flow  network  analysis 


The  preceding  equations  can  be  written  in  a  globally  linear¬ 
ized  form  with  Aj  being  the  globally  linearized  flow  coeffi¬ 
cients.  The  globally  linearized  loop  equations,  along  with 
the  inherently  linear  node  equations,  form  a  set  of  b  linear¬ 
ized  equations,  where  b  is  the  branch  number  (d  +  /).  In 
matrix  form,  this  can  be  represented  by: 

{e%^)lhi  =  B  0=h2,...,b)  (13) 

The  mass  flow  conservation  equations  are  linear,  but  the 
constitutive  equations  for  the  pressure  drops  are  only  linear 
for  low  Reynolds  number  (laminar  flow)  without  transitional 
losses.  In  such  cases,  Eq.  (13)  can  be  readily  solved  by 
the  matrix  method.  For  cases  involving  transitional  losses 
and/or  turbulent  flow,  a  system  of  non-linear  equations 
results,  which  can  be  solved  by  a  multi-variable  Newton’s 
method. 

The  global  linearization  method  used  to  solve  Eq.  (13)  is 
based  on  an  iterative  procedure.  The  computation  starts  from 
a  set  of  guessed  flow  rates,  which  determines  the  coefficient, 
Kj,  of  Eq.  (12),  and  solve  a  new  set  of  flow  rates  by  any  solver 
for  the  system  of  linear  equations.  An  iterative  procedure  is 


The  network  analysis  model  is  based  on  both  the  hydrau¬ 
lic  energy  principles  and  the  automatic  search  method.  The 
schematics  used  in  the  network  flow  code  are  presented  in 
Fig.  1 .  Channels  provide  a  passage  for  supplying  gases  to 
the  fuel  cells.  A  model  of  the  fuel  cell  flow  network  was  set¬ 
up.  A  multi-manifold  flow  scheme  was  used  to  ensure 
uniform  flow  through  the  different  gas  channels.  Gas 
stream  flows  into  the  manifold  and  three  cross-ducts  dis¬ 
tribute  the  flow  stream  into  the  trough.  Two  manifolds  were 
located  at  each  end  of  the  channels  and  connected  by  the 
cross-ducts. 

In  Fig.  1,  each  node  and  each  branch  are  numbered 
sequentially.  The  flow  directions  are  assumed  and  marked 
by  the  arrows.  If  the  flow  direction  inside  a  channel  is 
consistent  with  the  direction  assumed,  then  the  computed 
flow  rate  will  be  positive.  Otherwise,  the  sign  of  the  com¬ 
puted  flow  rate  will  be  negative.  For  a  manually  generated 
loop  equation,  the  loop  can  be  easily  identified  as  shown  in 
the  figure.  It  is  intuitive  to  choose  the  minimum  loop,  which 
includes  the  smallest  number  of  branches  and  the  indepen¬ 
dent  condition  for  the  loop  equation  will  be  automatically 
satisfied. 
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3.1.  Measurements  to  verify  the  network  analysis  scheme 

Experiments  were  run  to  measure  the  total  pressure  drop  in 
the  21  channel  flow  network  of  Fig.  1  in  order  to  check  the 
adequacy  of  the  various  empirical  loss  coefficient  data  used 
in  the  network  analysis  model.  An  experimental  apparatus 
was  designed,  built  and  instrumented  for  that  purpose  (Fig.  2). 

Referring  to  Fig.  2,  air  is  supplied  by  a  compressed  air 
cylinder  fitted  with  an  adjustable  pressure  regulator  and 
pressure  gauge.  The  volumetric  flow  rate  is  measured  using 
a  pair  of  rotameters  (Brooks  Instruments),  one  for  low  flow 
rates,  and  the  second  for  high  flow  rates.  A  three-way  valve 
allows  the  user  to  easily  switch  between  the  two  rotameters. 
Air  pressure  at  the  inlet  to  the  test  block  is  regulated  using  a 
test  pressure  gauge.  Flow  rate  and  back  pressure  are  con¬ 
trolled  by  a  pair  of  needle  valves,  one  at  the  inlet  to  the  test 
block,  and  the  other  at  the  exit.  The  pressure  drop  across  the 
test  block  is  measured  using  a  Dwyer  Microtector,  capable 
of  detecting  a  pressure  change  down  to  0.0002  in.  of  water 
(0.05  Pa). 

Experiments  were  run  at  regulated  pressures  of  193.1  kPa 
to  measure  the  total  pressure  drop  in  the  flow  network 
system.  In  the  flow  network  analysis  computation  using 
both  the  GT  code  and  SINDA/FFUINT,  the  total  pressure 
drop  from  the  network  inlet  to  its  outlet  was  also  calculated. 
Comparisons  between  measurements  and  calculations  are 
shown  in  Fig.  2. 

In  the  laminar  fully  developed  flow  regime,  the  friction 
coefficient  is  inversely  proportional  to  the  Reynolds  number. 
Therefore,  the  friction  loss  along  the  flow  path  varies 
linearly  with  the  increase  in  the  fluid  flow  rate.  However, 
the  transitional  loss  through  a  junction  is  expressed  as 
Cfpu2/2,  where  a  loss  coefficient  is  defined  as  Cf.  The 
pressure  drop  through  a  junction  increases  with  the  flow 
rate  in  a  square  power  form.  Hence  the  total  pressure  drop 


through  the  flow  network  varies  with  the  flow  rate  in  a 
quadratic  form. 

The  SINDA/FFUINT  analysis  does  not  include  the 
empirical  relations  for  the  transitional  losses  corresponding 
to  each  flow  junction  conditions.  Instead,  the  recommended 
constant  value  for  a  conjunction  was  used.  Therefore,  the 
predicted  SINDA/FFUINT  pressure  drop  curve  appears 
lower  and  more  linear  than  the  GT  code  prediction.  The 
experiments  show  a  strong  non-linear  effect  due  to  the 
quadratic  variation  of  transitional  losses.  The  simulated 
results  presented  in  Fig.  3  show  an  acceptable  agreement 
with  experimental  data;  the  calculated  error,  with  reference 
to  the  total  pressure  difference  along  the  manifold,  is  about 
4%. 

Another  experiment  was  performed  to  ascertain  if  laser 
Doppler  velocimetry  (FDV)  could  be  used  to  measure 
detailed  flow  velocity  distributions  in  the  various  channels. 
FDV  measurements  are  based  on  the  principle  that  two 
intersecting  laser  beams  produce  a  fringe  pattern  within  a 
very  small  target  area.  Pulses  of  light  reflected  off  particles 
passing  through  the  fringe  pattern  are  detected  by  a  photo¬ 
multiplier.  Signals  from  the  photomultiplier  are  sent  to  a 
signal  processing  system,  which  converts  the  pulse  rate  to  a 
velocity  value.  Given  a  sample  size  of  several  thousand 
particles  passing  through  the  target  area,  a  velocity  distribu¬ 
tion  can  be  produced  and  the  average  velocity  of  a  fluid  may 
be  determined. 

Experimental  error  in  the  FDV  measurements  is  largely 
due  to  the  many  factors  that  can  affect  the  measurements. 
Most  important  among  these  factors  is  the  placement  of  the 
target  beam.  Other  factors  include  small  fluctuations  in 
the  flow  rate,  electrical  noise,  reflectivity  of  the  surface  of 
the  channel,  density  of  the  seeding  particles,  and  the  presence 
of  tiny  bubbles  in  the  channel.  It  is  concluded  that  while  the 
FDV  measurements  are  useful  for  identifying  trends,  extra- 


Total  air  mass  flow  rate  (kg/s) 


Fig.  3.  Measured  vs.  predicted  pressure  drop  across  cathode  at  an  inlet  pressure  of  193.1  kPa. 
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Fig.  4.  Flow  distributions  through  channels  by  LDV  measurements,  along  with  calculations  from  SINDA/FLUINT,  GT  network  analysis  code,  and  the 
automatic  flow  network  generation  analysis  code. 


ordinary  measures  are  needed  to  assure  accuracy  of  the 
velocity  measurements  for  the  geometries  and  flow  condi¬ 
tions  of  interest  in  this  study.  The  care  and  attention  to  detail 
taken  in  this  study  resulted  in  velocity  distributions  consis¬ 
tent  with  prediction  of  the  codes  previously  validated  by  the 
pressure  drop  measurements. 

In  Fig.  1  the  channel  number  Cl  denotes  the  bottom 
channel,  and  the  channel  number  C21  indicates  the  top 
channel  of  the  schematic  flow  network.  Fig.  4  displays 
the  flow  distributions  through  the  21  channel  in  the  flow 
network  shown  in  Fig.  1  comparing  LDV  measurements,  the 
predicted  normalized  flow  rate  from  SINDA/FLUINT,  GT 
network  analysis  code,  and  automatic  network  equation 
generation  code. 

In  the  LDV  experiment,  the  pressures  at  the  inlet  manifold 
and  at  the  outlet  remain  constant  for  a  certain  flow  rate.  As 
exhibited  in  Fig.  4,  the  LDV  experimental  data  show  a 
certain  scatter  due  to  the  instabilities  of  the  decelerated 
flow.  The  simulated  results  present  an  acceptable  agreement 
with  experimental  data;  the  error,  which  is  calculated  with 
reference  to  the  channel  with  the  highest  deviation,  is  about 
4.7%.  Thus,  it  can  be  concluded  that  the  flow  in  the  inlet 
manifold  is  sufficiently  well  simulated  by  the  GT  code.  The 
result  by  automatic  network  equation  generation  agrees  well 
with  the  GT  code  as  well  as  the  LDV  experimental  data.  The 
SINDA/FLUINT  result  follows  the  flow  distribution  trend 
well  compared  with  the  LDV  measurement,  but  the  pre¬ 
dicted  distribution  is  more  uniform  than  the  actual  measure¬ 
ment.  This  error  is  mainly  caused  by  the  assumption  of  a 


constant  loss  coefficient,  which  is  difficult  to  modify  in  the 
SINDA/FLUINT  input  file  based  on  the  connection  and  flow 
conditions. 

In  the  analyzed  flow  network  shown  in  Fig.  1,  gas  stream 
flows  into  the  manifold  and  three  cross-ducts  distribute  the 
flow  stream  into  the  distribution  manifold.  Two  manifolds 
are  located  at  each  end  of  the  channels  and  connected  by 
the  cross-ducts.  The  cross-ducts  correspond  to  channels 
C4,  Cll,  and  Cl 8.  From  Fig.  4,  the  results  of  the  LDV 
measurement  show  peaks  in  those  three  channels.  Channel 
C4  is  closest  to  the  inlet  and  outlet  side  and  has  highest  flow 
rate.  The  flow  rate  trend  decreases  away  from  the  inlet  and 
outlet  side,  which  indicates  flow  favoring  the  shortest  short 
path. 

The  design  of  the  two-manifold  connected  by  cross-ducts 
to  distribute  inlet  or  exit  gas  provides  a  way  to  adjust  the 
peak  flow  position  and  also  the  flow  distribution  pattern.  For 
achieving  an  optimum  flow  distribution,  the  flow  network 
design  is  essential.  With  the  aid  of  the  flow  distribution 
analysis  through  computation,  the  selection  of  the  manifold 
size  and  network  configuration  for  realizing  a  desired  flow 
distribution  become  easy.  In  a  flow  network  design,  better 
flow  uniformity  may  be  achieved  if  the  ratio  of  the  pressure 
drop  across  the  channels  to  the  pressure  drop  in  the  plenum 
and  other  connecting  ducts  is  maximized,  without  significant 
increases  in  the  overall  pressure  drop.  One  design  approach 
leading  to  a  uniform  flow  distribution  is  to  minimize  the 
pressure  drop  in  the  plenum  without  significantly  increasing 
the  pressure  drop  in  the  channels. 
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4.  Conclusions 

The  methodology  used  in  the  GT  code  for  flow  network 
analysis  has  been  presented.  The  code  includes  transitional 
losses,  based  on  the  loss  coefficients’  data  in  Idelchik’s  flow 
handbook.  The  code  can  be  used  to  describe  the  flow  distribu¬ 
tion  through  a  complex  network.  The  code  provides  added 
flexibility  in  modeling  such  networks,  compared  to  SINDA/ 
FLUINT,  the  commercially  available  flow  network  analysis 
code.  An  automatic  loop  identification  and  network  equation 
generation  code  has  been  developed  and  shown  to  be  con¬ 
sistent  with  the  manually  generated  equations  in  the  GT  code. 
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